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THE  BRAODENING  OF  SPECTRAL  LINES  BY  AUTOIONIZATION, 

RADIATIVE  TRANSITIONS,  AND  COLLISIONS 

I.  INTRODUCTION 

The  shape  of  atomic  spectral  lines  is  seldom  determined  exclusively  by 
the  spontaneous  radiative  or  autoionization  processes.  In  addition  to  the 
Doppler  effect,  spectral  lines  can  be  appreciably  broadened  as  a  result  of 
collisions  with  the  surrounding  particles,  classical  electromagnetic  fields, 
and  radiation  fields.  Although  natural  broadening  is  neglected  in  most 
theories  of  collision  or  pressure  broadening,  it  is  desirable  to  have  a  com¬ 
prehensive  theoretical  framework  in  which  all  broadening  mechanisms  can  be 
treated  on  an  equal  footing.  Such  a  theory  could  be  applied  to  calculate  the 
shapes  of  satellite  lines  which  are  associated  with  resonance  lines  of  multiply- 
charged  ions  in  the  x-ray  emission  spectra  of  high-temperature  plasmas.  Since 

the  satellite  lines  are  prominent  features  over  a  wide  range  of  electron  den- 
10  24  -3 

sities  (10  -  10  cm  ) ,  it  is  desirable  to  take  into  account  their  natural 

width  due  to  spontaneous  radiative  decay  and  autoionization.  Another  appli¬ 
cation  of  a  more  general  theory  would  be  the  simultaneous  treatment  of  colli¬ 
sion  and  radiation  broadening. 

Most  treatments  of  collision  or  pressure  broadening  are  based  on  the 
evaluation  of  the  autocorrelation  function  of  the  dipole-moment  operator.  The 

modem  quantum  theory  of  spectral  line  broadening  by  plasmas  has  been  developed 

1-3  4 

from  this  point  of  view  by  Baranger  and  independently  by  Kolb  and  Griem  . 

Collision  broadening  by  neutral  gases  has  been  treated  by  Fano”*  using  a  varia¬ 
tion  of  the  autocorrelation-function  approach  which  is  directly  applicable  to 

other  broadening  mechanisms.  The  key  quantity  in  Fano's  treatment  is  the 

£ 

frequency-dependent  relaxation  operator  introduced  by  Zwanzig  .  This  operator 

obeys  a  Lippman-Schwinger  type  equation  which  is  defined  in  the  Liouvi lie-space 
Manuscript  submitted  September  2,  1980. 
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representation,  in  which  density  matrices  are  described  by  vectors.  In  the 

b 

present  investigation,  we  derive  the  general  form  of  the  overlapping  line-shape 
1-4 

formula  by  evaluating  the  Liouvllle-space  expression  for  the  spontaneous 

electric-dipole  emission  rate  between  states  of  the  complete  quantum  mechanical 

system,  which  consists  of  the  atomic  system,  the  surrounding  particles,  and  the 

radiation  field.  The  key  quantity  in  this  treatment  is  the  Liouville-space 

scattering  operator,  employed  in  the  investigations  of  Ben-Reuven  and 
8~10 

co-workers  ,  which  is  defined  in  terms  of  the  Zwanzig  relaxation  operator. 

The  determination  of  the  spectral  line  shape  is  reduced  to  the  problem  of 
evaluating  the  matrix  elements  of  the  relaxation  operator,  whose  real  and 
imaginary  parts  correspond  to  the  line  shifts  and  widths,  respectively.  In 
this  investigation,  the  explicit  expressions  for  the  widths  are  obtained  in  the 
isolated-line  approximation  by  evaluating  the  lowest-order  nonvanishing  contri¬ 
bution  to  the  diagonal  matrix  elements  of  the  relaxation  operator.  We  show 
that  the  total  width  is  the  sum  of  the  rates  for  all  inelastic  transitions  from 
the  initial  and  final  atomic  states  comprising  the  spectral  line  and  a  term 
involving  the  square  of  the  difference  between  the  elastic  scattering  amplitudes. 
In  a  treatment  of  the  isolated  line  shape  problem  based  on  the  autocorrelation 
function,  Lambropoulos^  obtained  only  the  inelastic  contributions  and  did  not 
allow  for  autoionization  processes. 

The  remainder  of  this  paper  is  organized  as  follows:  In  section  II,  we 
review  the  general  theory  of  transition  probabilities  within  the  framework  of 
the  density  matrix  description  of  atomic  states.  In  section  III,  the  spontaneous 
emission  rate  for  electric-dipole  radiation  is  evaluated  and  shown  to  have  the 
same  form  as  the  general  formula  for  overlapping  lines.  In  section  IV,  explicit 
expressions  are  obtained  for  the  widths  associated  with  various  broadening 
mechanisms  in  the  isolated-line  approximation.  The  variations  of  these  widths 
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with  the  atomic  states  and  with  the  properties  of  the  surrounding  particles  and 
radiation  field  are  discussed  in  section  V. 

II.  DENSITY  MATRIX  DESCRIPTION 

A  very  general  theory  of  transition  probabilities  can  be  developed  within 
the  framework  of  the  density  matrix  description  of  states  in  quantum  mechanics 
In  order  to  describe  radiative  transitions  of  an  atomic  system  which  is  inter¬ 
acting  with  the  surrounding  particles  through  collisions,  it  is  necessary  to 
first  consider  the  total  Hamiltonian 

H  =  HS  +  HP  +  HR  +  VSP  +  VSR  +  VPR  ,  (1) 

S  P 

where  H  is  the  Hamiltonian  for  the  isolated  atomic  system,  H  is  the 

Hamiltonian  for  the  perturbing  particles,  and  Hv  is  the  Hamiltonian  for  the 

free  radiation  field.  The  interaction  between  the  atomic  system  and  the 

SP 

surrounding  particles  is  represented  by  V  ,  while  the  interaction  between 

SR 

the  atomic  system  and  the  radiation  field  is  denoted  by  V  .  The  interaction 
PR 

V  , which  describes  radiation  processes  involving  the  perturbing  particles,  is 
usually  neglected  in  the  theory  of  spectral  line  shapes.  However,  a  treatment 

of  perturber  radiation  interfering  with  line  radiation  has  been  given  by 

D  12 
Burgess 

Eigenstates  of  H^,  HP,  and  HR  will  be  denoted  by  |a\  |b>,  ...,  |p>, 

|q>,  ....  and  |n\  |m>  ...,  respectively.  The  index  p  represents  the  quantum 
states  of  all  the  surrounding  particles,  and  the  index  n  is  understood  to 
represent  a  set  of  photon  occupation  numbers,  one  for  each  mode  of  the  radi¬ 
ation  field.  The  notation  can  be  further  compressed  by  introducing  the  single 
index  a  to  represent  a,  p,  n;  and  the  index  £  for  b,  q,  m,  etc.  We  shall  assume 
that  the  zero-order  Hamiltonian  eigenvalue  problem 

H°  j  a)  ■  #i  u  |  a)  >  (2) 

i  a  ' 


where 


can  be  solved.  The  solution  of  (2)  pan  be  most  easily  obtained  If  correlations 

between  the  perturbing  particles  are  neglected.  In  the  theory  of  spectral  line 

broadening  by  plasmas4"  ,  |a>  and  |b>  are  often  taken  to  be  eigenstates  of  the 

atomic  system  In  the  presence  of  the  electric-microfield  produced  by  the  plasma 

ions,  which  is  usually  assumed  to  be  time-dependent  on  the  time  scale  of  the 

radiative  transition.  In  order  to  describe  autoionization,  part  of  the  electro- 

5 

static  interaction  between  the  atomic  electrons  must  be  removed  from  H  and 
included  in  equation  (6)  below. 

A  general  eigenstate  of  H°  can  be  represented  by  the  density  operator 


S  P  R 
P±  =  P  x  p  x  p 


(4) 


for  the  initial  state.  The  separable  form  implies  that  initially  there  are  no 

correlations  between  the  atomic  system,  the  perturbing  particles,  and  the  radi- 

S  P  R 

at ion  field.  The  matrices  corresponding  to  p  ,  p  ,  and  p  are  assumed  to  be 
diagonal  in  the  ^presentations  a,  p,  and  n,  respectively.  For  example,  the 

S 

atomic  density  matrix  p  can  then  be  expanded  in  the  form 

PS  “  £  P„  |a>  <  a|  .  (5) 

a  a 

In  order  to  describe  the  transitions  between  the  eigenstates  of  H°  which 
are  induced  by  the  interaction 


V 


VSP  +  VSR  +  V**  , 


(6) 


it  is  necessary  to  investigate  the  time  evolution  of  the  density  operator 
p(t)  for  the  complete  interacting  system.  After  Introducing  the  Liouville- 
space  operator 

(H  I*  -  I  H*), 

the  Schrodinger  equation  can  be  written  in  the  form'* 


(7) 


(8) 


—  p(t)  “  -  i 


dt 


p(t) . 


The  matrix  elements  of  equation  (8)  are  given  by 


pa  b  (t) 


ab,cd  pcd(t) 


where 


2,  .  ,  =  ft”1  (H  6,  ,  -  <5  H*  )  . 

ab,  cd  acbd  ac  bd 


(9) 


(10) 


The  asterisk  denotes  complex  conjugation.  The  solution  to  equation  (8)  can 
be  expressed  formally  by 

P  (t)  =  V  (t,tQ)  p(tQ)  ,  (11) 

where  V  (t,t  )  is  the  Liouville  time-evolution  operator  defined  by 


y  <t.t  ) 


=  exp 


[-it 


(12) 


In  order  to  define  the  transition  probability,  it  is  necessary  to  intro¬ 
duce  the  final  state  projection  operator  P^.  The  spontaneous  emission  of  a 
single  photon  with  wave  number  k  and  polarization  X  can  be  described  by  the 
final  state  projection  operator 

Pf=/CU><b|x  |l^A><lj>JxIB,  (13) 

b 

where  Ig  is  the  unit  operator  in  all  other  degrees  of  freedom  which  will  be 
treated  as  bath  states,  to  be  defined  below.  The  probability  per  unit  time 
for  the  transition  i  f  can  then  be  expressed  by^ 


A(i  -*•£)  -  lim  lim  Tr 

00  £  OO 

o 


<t.to)  pt  <t0)| 


(14) 
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^•n  ^ 


which  is  the  Liouville-space  analogue  of  the  definition  given  by  Goldberger 
and  Watson^. 

Equation  (14)  provides  a  general  description  of  a  collision  process  in 

which  the  asymptotic  states  are  eigenstates  of  the  Hamiltonian  H°  for  the  non- 

9 

interacting  systems.  It  has  been  pointed  out  that  the  same  description  can 
be  applied  to  relaxation  processes  which  involve  an  initial  stationary  state 
(e.g.  a  state  with  only  diagonal  non-zero  density  matrix  elements). 

The  asymptotic  time  limits  can  be  taken  to  obtain  the  steady-state 
transition  rate^ 

A(i  -f)  =  -  i  Tr  |  Pf  T(+  10)  Pi|  .  (15) 

The  Liouville  operator  ^^(z)  is  the  tetradic-space  analogue  of  the  T  matrix 
occurring  in  ordinary  scattering  theory.  The  Liouville-space  analogue  of  the 
Lippman-Schwinger  equation  can  be  written  in  the  form 

^(z)  *  'lf+'V/7£(z)  V  ,  (16) 

where ^3L(z)  is  the  Green's  operator 

faz)  -  (z  -  &)'1  .  (17) 

The  Liouville  operator  has  been  divided  according  to 

<£.  »,<L0  AT,  (is) 

where  and  V  are  the  Liouville  operators  corresponding  to  H°  and  V, 

respectively.  The  limit  z  -*•  i  0  of  the  complex  variable  z  in  equation  (15) 
corresponds  to  reaching  the  origin  from  above  the  real  axis.  The  initial  state 
density  operator  and  the  final  state  projection  operator  in  equation  (15) 

must  correspond  to  stationary  states  of<J-°,  i.e. 


X.°  »i-  *£- 


(19) 


In  the  density  matrix  formulation  of  relaxation  processes**^'10  ,  the 
transition  amplitude  is  expressed  in  terms  of  a  scattering  operator*X  r  (z), 
which  is  defined  on  a  subspace  of  states  incorporating  only  the  interaction 

—V 

between  the  uncorrelated  atomic  system  and  the  relevant  mode  k  \  of  the 
radiation  field.  The  scattering  operator  *"Xr(z)  is  then  explicitly  expanded 
in  a  power  series  of  the  relevant  interaction/VPr ■  The  remaining  interactions 
are  introduced  as  relaxation  effects  by  a  self-energy  operator.  The  distinction 
between  the  "free"  states  of  the  relevant  degrees  of  freedom  and  the  states  in 
which  they  are  correlated  to  the  remaining  (bath)  degrees  of  freedom  can  be 
formally  achieved  by  using  the  Liouville-space  projection  operators  TP  and 
GUi-TP  which  were  introduced  by  Zwanzig^.  These  projection  operators 
play  a  role  similar  to  that  of  the  dyadic  projection  operators  introduced  in 
the  theory  of  nuclear  reactions  by  Feshbach^.  Their  precise  definition  will 
be  given  below  later. 

The  relevant  part  of  the  interaction  can  be  defined  by 

If'  -  tP'Vf'lR.  (20, 

The  relevant  scattering  operator  *X  r(z)  obeys  the  Lippmann-Schwinger  type 
equation 

where 


[z  - -Pt°?P- TP 'l  • 

The  self-energy  operator  0^(z)  is  given  by 

*<«>-  vau-aia  r1 


(22) 


(23) 


Equations  (22)  and  (23)  are  analogous  to  the  results  obtained  for  the  ordinary 

(dyadic)  Greens  operator  by  Mower^.  The  transition  rate  can  be  evaluated 
in  terms  of  f(+  i  0  )  as  follows^. 


_7_ 


(24) 


A(i  +  f)  -  -  i  Tr  |pf  t(+  i  0)  p± 

The  Zwanzig  projection  operator  can  be  defined  in  the  double-space  notation 
1-3 

introduced  by  Baranger  as  follows: 

IP  =  |pB»«IB|,  (25) 

where  B  denotes  the  bath  degrees  of  freedom.  The  effect  of^^Ron  a  Liouville 
operator  is  to  average  it  over  these  degrees  of  freedom.  The  double-space 
vectors  in  the  subspace  spanned  by  IP  will  be  denoted  by  |  ab,  n  m  >)>  ,  where 
|  a  b  )  )  represents  the  operator  |  a  )>  <  a  |  in  the  space  of  eigenstates  of 

g 

H  and  |  n  m  >  )  represents  |  n  >  <  m  |  in  the  space  of  relevant  photon  occupation 

j- 

number  states.  The  initial  state  density  operator  and  the  final  state 

projection  operator  Pj  ,  which  are  defined  on  the  relevant  degrees  of  freedom, 
can  be  expanded  in  the  form 

=  L  p  |  a  a,  0  0  >>  (26) 

a  a 

and 

Pj  =  ^  |  bb,  1  1  »  .  (27) 

The  single-photon  spontaneous  emission  rate  can  now  be  expressed  in  terms 
of  the  matrix  elements  of  r(+  i  Oj  as  follows: 

A  (i  f )  =  -  i  <  <  Pf  |fXr(+iO)|  p±>> 

=  -i  EZ  <<  bb,  11  |  r  (+  i  0 )  |a  a ,  00»p  .  (28) 

a  b  a 

The  summations  over  a  and  b  are  to  be  taken  over  all  bound  and  continuum  states 
of  the  atomic  system.  In  practice,  these  summations  are  restricted  to  include 
only  a  few  discrete  states. 
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III.  SPONTANEOUS  EMISSION  OF  ELECTRIC-DIPOLE  RADIATION 


In  this  section, we  evaluate  the  spontaneous  emission  rate  for  single¬ 
photon  electric-dipole  radiation.  The  objective  of  this  evaluation  is  to 
demonstrate  that  the  formula  for  overlapping  lines^  ^ ,  which  is  usually 
obtained  by  starting  with  the  autocorrelation  function  of  the  electric-dipole 
moment  operator,  can  be  derived  from  the  general  equantum  mechanical  expression 
for  the  transition  rate  ir;  the  Liouville-space  representation  of  density 

matrices.  Alternative  density  matrix  formulations  of  the  line  shape  problem 

17  18 

have  been  given  by  Cooper  and  Davis  .  In  these  formulations  the  time 
dependence  of  the  density  matrix  is  treated  directly. 

Since  the  second  term  in  equation  (21)  gives  the  only  nonvanishing  con¬ 
tribution  to  the  scattering  operator  “ji  r(z),  the  rate  for  the  spontaneous 
emission  of  dipole  radiation  is  given  by 


A  (w) 
r 


„  _  j  £  £ 

a  b 


bb 


.  ii  nr&r  (+io>  nn  aa.oo 


-i  £  £  0  £  £  £  £ 

a  b  a  cn  d  n ,  e  n  f  n  . 

c  a  e  t 

bb,  11  |  I  cd>  nc  nd 


kX 


e  f  ,  n  nf 


nc  nd  |Jlr<+  1  °> 

\$1‘  »e  nf  I'XT kt  I  a  *•  °  °^> 


(») 
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where  all  summations  over  the  relevant  degrees  of  freedom  are  indicated 
explicitly  in  the  second  version. 

The.  matrix  elements  of  the  Liouville  interaction  operator  XT 

fx 

can  be  evaluated  using  the  definition  which  follows  from  equation  (7).  The 
matrix  elements  of  the  dyadic  interaction  operator  V  are  evaluated 


using  the  expansion 


vix  I  a> n  >’ 


Vb  a  (k) 


+  Vba< 


t)  ^n  j  b,  n  -  1^ 

-£)^n  +  L  |  b,  n  +  l^J  , 


where 


Vba(k)  ' 


■-  S  f  v'j 


We  then  obtain  the  result 


,  ne  |  *  01))  ’ 

fi‘1|  %  (-k*>  5n  1  6  f  «  0  1 

-  6  6  n  V*  (-£)  6  i  ?  . 

e,  a  ng,  0  fa  nf  l  ) 

We  obtain  a  similar  expression  for  <<  bb,  11  |  V?  i  cd,  nc  nd  >>  , 

except  for  the  presence  of  terms  corresponding  to  n  or  n.  ■  2.  These  terms 

c  a 

can  be  neglected  because  they  do  not  contribute  to  the  dominant  spontaneous 
emission  process. 

The  frequency-dependent  transition  rate  Ar(w)  given  by  equation  (29)  can 
now  be  written  as  a  sum  of  four  terms.  These  terms  may  be  combined  in  pairs 
by  utilizing  the  Hermitian  property 

*ab  <k'>  ■  V.  <k) 
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and  the  Liouville  conjugation  symmetry  relationship 


^ab’  cd  (z) 


‘  ^ba’  dcM> 


which  is  a  consequence  of  the  microscopic  reversibility  of  the  time  evolution 
operator  under  time  reversal. 

The  frequency-dependent  transition  rate  can  now  be  written  in  the  form 


Ar(u>) 


=  _  2_  £  £  £  £  v 

*2  im  a  b  c  d  db 


(w)  V  k)  p 


bd,  ac  ac 


=  -  ^  Im  Tr  [  *(k) /£ 


(go)  p  V  (k) 


where  the  reduced  Green's  operator  (go)  is  defined  by 

d  ,  ac  (w)  =  <^b  d  *  01  l^r(+io)  I  ac’° 

(36) 

-  ^<^bd,  1  0  lA  r<+  1  °>  lac*°  1^- 

The  second  term,  in  which  the  photon  occupation  numbers  are  interchanged  in  the 
bra,  does  not  give  a  sharp  resonant  contribution  to  the  spectral  line 
shape  and  will  therefore  be  omitted. 

Using  the  relationship 

?  3L°  HP  I  a  b,  n  m  *  £wa  -  u>b  +  (n-m)  go]  |  ab,  nm^  (37) 

and  neglecting  the  relaxation  operator  R  ,  the  reduced  tetradic  Greens  operator 
generally  defined  by  equation  (36)  has  matrix  elements  given  by 


b  d ,  a  c 


(go)  ■  lim  (go  -  go  +  go  +  i  n)  6.  _  6 


b  a  d  c 


Consequently,  the  frequency-dependent  transition  rate  in  the  absence  of  relaxation 
effects  simplifies  to 


11- 


(39) 


A  r(w) 


2jr 

2 


-ft 


E 

a 


r 

b 


a  b 


(k) 


6  (u>-w  + 

a 


G)b) 


which  is  the  familiar  Fermi  Golden  Rule  formula  for  spontaneous  emission  from 
the  levels  a  to  the  levels  b. 

In  order  to  obtain  the  differential  transition  rate  for  the  emission  of 

a  photon  per  unit  solid  angle  and  angular  frequency  intervals,  equation  (35) 

2  3 

must  be  multiplied  by  the  density  of  final  states  factor  to  /(2rrc)  .  This  factor 
would  have  been  introduced  automatically  if  the  final  state  projection  operator 
(27)  had  been  more  carefully  defined  to  take  into  account  the  continuous  distri¬ 
bution  of  modes  of  the  radiation  field.  Using  the  electric-dipole  approximation 


V  (k) 
b  a 


=  -  i^TTlT^  ^b  |  Dx  |  , 


(40) 


where  denotes  the  component  of  the  total  atomic  dipole-moment  operator  in  the 

direction  of  polarization,  the  differential  photon  emission  rate  for  a  given 

polarization  X  is  obtained  in  the  form 

3 


A  (w,Q)  ■ 


2ir  ft 


“3  L  A  (u»°  ) 


(41) 


The  line  shape  function  L.  (w,ft),  which  depends  on  the  angles  0  of  photoemission. 


is  given  by 


L,  (w,n  )  =  -—  Im  Tr 

A  TT 


k&<“>  <">»] 


(42) 


Usually,  one  is  Interested  only  in  the  total  spontaneous  emission  rate  per 
unit  frequency  interval 

/•  ,  /.  ..3 

(43) 


A  (to) 
r 


E 

X 


where 


L(u>)  ”  -  —  Im  Tr 

IT 


f  dn  a  A  (u>,  n)  -  ^  l(u>), 

J  "ft  c  J 

[  eT  •  p  •  D  ]  . 


(44) 
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Equation  (44)  has  the  same  form  as  the  results  derived  by  Baranger  ,  by 

Kolb  and  Griem^,  and  by  Fano^.  Note,  however,  that  equation  (44)  is  more 
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general  than  the  formula  usually  employed  for  overlapping  lines  in  that  no 
restriction  has  as  yet  been  imposed  on  the  subspace  of  the  atomic  states  over 
which  the  trace  is  to  be  taken.  Consequently,  equation  (44)  provides  a 
general  framework  for  the  description  of  the  entire  emission  spectrum  due  to 
all  bound-bound  transitions,  including  transitions  involving  autoionizing  levels. 


IV,  The  Isolated  Line  Approximation 

If  only  the  diagonal  matrix  elements  a  b  ,  ab  (uj)  of  the  reduced 
Green's  operator  U  (w)  are  included  in  the  line  shape  formula  (44) ,  we  obtain 


L (to)  =  -  —  Im  ^  ? 

u  a  b 


Dab 


W-0)  +  w  -19  h. 


ab  *  a  b 


(u>) 


(45) 


where 


'ft  a  h  .  h  <«)  - 

3D,  3D 


^ab,  0  i  |  T'ft  (+  i  0)  1P\  ab,  o)).  (46) 

From  the  definition  of  the  tetradic  relaxation  operator  ft  (z)  given  by  equation 
(23),  it  follows  that 

*X+(z)  -  'ft  (z*)  ,  (47) 

where  the  Hermitian  conjugation  of  ft  is  defined  by 

<Ti+)„e,  ye  -n* .  (48) 
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Consequently,  the  procedure  of  Goldberger  and  Watson  can  be  followed  to 


obtain 


(49) 


^“o  1R.  (X  ±  i  n)  -  A(x)  +  i  r(x)  I  2, 

where  A(x)  and  r(x)  are  Hermitian  operators, 
expressed  in  the  form 


ab,ab 


(to)  can  therefore  be 


^  ab  .a  b(u) 


A  .  (w)  -  i  T  ,  (a.)  2, 

a  b  ,a  b  a  b,ab 


(50) 


where  A  ,  ,  (to)  gives  the  shift  associated  with  the  line  a  b  and 

3l  b  j  fl  b 

r  3  b  a  ^  (<*>)  corresponds  to  the  full-width  at  half  maximum. 

It  Is  clear  that  equation  (45)  corresponds  to  an  incoherent  superposition 

of  the  Lorentzian  profiles  associated  with  each  atomic  line  a-*-b.  The  isolated 

line  approximation  is  expected  to  be  valid  provided  that  the  separations  of  the 

lines  are  large  compared  with  their  widths.  In  this  section,  we  evaluate  the 

lowest-order  nonvanishing  contributions  to  the  widths  arising  from  autoionization, 

spontaneous  and  induced  radiative  transitions,  and  collisions. 

The  width  and  shift  operators  T  (x)  and  A(x)  can  be  expressed  in  terms  of 

the  relaxation  operator  ^  (x  +  i  q)  by  following  the  procedure  used  by 

14 

Goldberger  and  Watson  .  From  the  definition  (23)  and  the  identity 

(z-Ql1!-(X  )_10 -  (z  -  Qi^OO"1  <X7l(z),  (51) 


which  is  analogous  to  the  dyadic  relationship  used  by  Mower1**,  we  obtain  the 
results 

T(x)  ■  2tt  lim  \tl  +  O'  +  i  n)  Qi  6  (  x  -  <3l  0  (fc.  )  QlT?  (x  +  i  n)1 

n  o  L  J 


(52) 


r  (X)  -  Tie  V  -  — 

2ir 


I 

/ 


and 

r(x»)  dx  1 
x’  -  X 


(53) 


where  P  denotes  the  Cauchy  principal  value.  Similar  results  have  also  been 
derived  by  Lambropoulos11. 
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The  lowest-order  nonvanishing  contribution  to  the  width  operator  TCx), 
which  is  obtained  from  the  approximation  of  ft  (z)  by  \f  ,  is 

r(z)  =  2  tt  1/  Q.  6(x  -  Q  ° 


(54) 


Using  the  compressed  notation  introduced  in  section  II  for  the  eigenstates  of 

'-f  o 

ot-  ,  the  general  matrix  elements  of  T(x)  can  be  expressed  in  the  form 


|„  6))  '  ?  5,.. 


<S(x-a)a, ,  +  w., ,) 


8  |VI*  -  IV*  |  a"  6 

2’v 


’  3 "  |  V  I*  -  I  V*  | 


2tt 

fi2  a" 


"  V  »•  •  6q  *  6(x  -  w  "  +  a»c) 

a  a  a  a  B  3  a  & 


+  2  tt 


4’  E 


ft2  6 


,,  V|  6"  6aa’  6(x"  \  +  V* 


(55) 


H  Vaa*  Vg8’  6(x  -V  +  V 

Tl 


V.a-  '•Ig  5  (*-“<,  +  V)  • 


where  we  have  made  the  approximation  Gc  =  1,  discussed  by  Ben-Reuven 
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To  obtain  the  line  width,  we  must  evaluate  the  expression 

r  .  .  M 

a  b  ,  a  b 


(^ab,  0  1  I  *2P  T  (+  i  0 )  IP  |  ab,  0  l\  , 


(56) 
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It  is  now  necessary  to  introduce  the  bath  degrees  of  freedom  n  and  p,  where  n 
represents  the  set  of  all  photon  occupation  numbers  and  p  is  the  quantum  state 
of  the  perturbing  particles.  The  frequency-dependent  width  can  then  be  expressed 
in  the  form 


r  ,  .  (uj)  - 

ab  ,  a  b 


Z  Z  Z  oi» n  n 

n  n'  p  p' 


pp  I  r(o)| 


(57) 


ab,  0  1,  n'n',  p'  p 


It  is  necessary  to  make  the  somewhat  artificial  distinction  between  the  relevant 
and  bath  photon  modes  in  order  to  obtain  the  frequency-dependence,  which  occurs 
in  the  delta  functions  in  equation  (55).  In  subsequent  equations,  reference  to 
the  relevant  photon  mode  numbers  will  be  omitted  whenever  it  becomes  unnecessary. 
After  substituting  equation  (55)  into  equation  (57),  the  lowest-order  non¬ 
vanishing  contribution  to  the  Isolated  line  width  can  be  obtained  in  the  form 


ab 


a  b  (<»)  -  2  2  2  ES  Kanp  |  V  |  a"  n"  p\ 

’  -ft  a»^a  n  n "  p  p M  ' 


x  6  (w  -  u  „  +  <•>.-«.,  +  «  -«,.+»)  p  p 
a  o  p  p  n  n  np 

^  ^  Z  Z  Z  Z  |  <(b  n  p  |V|  b"  n"  p"\ 

“ft  b"^b  n  n"  p  p"  / 


P  P 


x  6  (a)  -  o>a  +  -  <op  +  ay,  -  u)n  +  a.n„)  pn  pp 

+  ^  ^  I  ^JP  |v|  a  -  ^bp|v|bp^> 


(58) 


P  P 


x  6  (u  -  u>a  +  Wjj  -  u  .  +  «o  )  ppi 
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The  first  two  terms  describe  all  inelastic  transitions  from  the  initial  and 
final  states  of  the  radiating  atomic  system.  Note  that  the  states  of  the  radi¬ 
ation  field  are  not  involved  in  the  last  term,  which  describes  only  elastic 
collisions.  This  is  a  consequence  of  the  assumption  that  photons  are  emitted 
only  by  the  atomic  system  and  the  neglect  of  correlations  between  the  atomic 
system  and  the  perturbers.  The  radiation  emitted  by  the  perturbers  can  be 
described  by  the  use  of  higher-order  perturbation  theory. 

A.  Autoionization 

g 

Autoionization  results  form  the  electrostatic  interaction  V  between 
a  discrete  state  and  the  continuum  states  of  the  atomic  system  having  the  same 
energy.  The  definition  of  the  zero-order  atomic  states  can  be  made  in  an 
unambiguous  manner  by  using  the  dyadic  projection  operators  introduced  by 
Feshbach^"*  or  by  the  alternative  procedure  introduced  by  Fano2^.  Autoioni- 
zation  can  be  included  in  the  present  theory  by  transferring  V  from  H°  to  V, 
as  discussed  in  Section  II.  The  width  arising  from  the  possibility  of  auto¬ 
ionization  from  both  the  initial  and  final  states  of  the  radiating  atomic 
system  is  obtained  from  the  inelastic  terms  in  equation  (58)  and  can  be 
written  as 

rAab,ab  (w)  =  ^  I  <^a  I  VS  |  c^>  |2  6(«-uc  +  o>b) 

(59) 

+  ^  I  ^  b  j  VS  |  c^)>  | 2  6  (to  -  wa  +  wc) , 

where  the  sum  over  c  ranges  over  all  continuum  states  which  are  allowed  by  the 

energy  and  symmetry  selection  rules.  If  we  now  make  the  approximation 
w  ■  w  -  in  which  the  frequency-dependence  is  neglected,  we  obtain  the 
result 
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(60) 


which  is  the  sum  of  the  autoionization  rates  from  the  initial  and  final  atomic 
levels.  In  its  most  general  form,  equation  (60)  gives  the  autoionization  con¬ 
tribution  to  the  line  width  for  a  radiative  transition  between  two  autoionizing 
levels  a  and  b. 

It  is  well-known  that  the  mixing  between  the  discrete  and  continuum  states 

20 

produces  an  asymmetric  profile  .  However,  to  the  best  of  the  authors  knowledge, 

no  theories  have  been  developed  which  incorporate  the  effects  of  radiative 

transitions  and  collisions  into  the  asymmetric  line  shape  function.  In  the  limit 

20 

where  the  line  profile  parameter  q  introduced  by  Fano  becomes  large,  the  line 
shape  function  approaches  a  Lorentzian  function,  which  is  generalized  in  the 
present  investigation  in  order  to  treat  all  broadening  mechanisms  on  an  equal 
footing. 

B,  Radiative  Transitions 

The  width  arising  from  all  radiative  transitions  from  the  states  a 
and  b  which  is  obtained  from  the  inelastic  terms  in  equation  (58)  is  given  by 


x  6  (w  -  w  +  w.  i,  -  oo  +  w  ,,)  p_ 
ao  n  n  n 
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We  now  introduce  the  expansion  of  the  interaction  V  in  terms  of  modes  kX  of 

->■ 

the  radiation  field  and  replace  the  summation  over  k  in  this  expansion  by  inte¬ 
grations  over  the  frequency  and  the  solid  angle  52^.  In  the  electric-dipole 
approximation,  the  total  radiative  contribution  to  the  width  is  given  by 

rRab,  ab  <“>  '  — ,  ^  ff  “k  d"k  d0k 
2  uh  c  X  J  J 

(ua  ~  wa»)2  (" 5  (w  ~  wa"  +  ^  “  wk) 


+  <^1^  (ojk,  f2k)  /  6  (w  -  a>a  +  +  a)k) 

+  <CX  (V  o,k^>  6  (to  -  u>a  +  %n  -  V  j 


where  <  n^  (u>k,  >  denotes  the  average  number  of  photons  per  unit  frequency, 

per  unit  solid  angle,  which  have  polarization  X. 

The  contribution  from  all  spontaneous  radiative  transitions  can  be  obtained 
from  the  terms  which  are  independent  of  the  average  number  of  photons  by  per¬ 
forming  the  integration  over  &k  and  the  sum  over  X.  If  we  now  make  the  approxi¬ 
mation  a)  *  u)  -w.  ,  we  obtain  the  result 
a  b 


ab’  8  b  3*  c3  a"  +a  a  a 


+  3-fic3  b"  +  b  (a)b  ‘  V >' 


<**«  ~  “  i.)  <  a  D  a' 


<(b  I  D  |  b ”>  | 


e(w  -  u>  ,,) 

ct  a 


0(ub  -  V)! 


where 


l1  K  >  1  ) 
(0  x  <  1  ) 


which  is  the  sum  of  all  spontaneous  emission  rates  out  of  the  states  a  and  b 
comprising  the  spectral  line. 

The  contribution  from  all  Induced  radiative  transitions  is  conventionally 
expressed  in  terms  of  the  specific  intensity  I^(to,  )  which  is  related  to  the 
average  number  of  photons  by 

3  2 

<  nA  (a.,/!)  >  -  Ix  <to,A).  (65) 

u 

If  we  again  make  the  approximation  to  “  to^  -  we  obtain  the  result 


b  ,  a  b  72 


tic  X 


I  !“  4"/>  I  [  VV“a  'JV0<V“a  >  0  <V’"“a)] 


♦  5  1 4 1.  i2  [x 

b»^b  1  >  /  [ 


1 x<<VV/lk)G  C^u-tob) 


which  gives  the  contribution  to  the  width  from  all  stimulated  radiative  emission 
and  absorption  processes  out  of  the  states  a  and  b.  Results  similar  to  those 
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given  by  equations  (63)  and  (66)  have  also  been  obtained  by  Lambropoulos 
The  starting  point  in  his  work  is  the  autocorrelation  function. 


C.  Collisions 


The  contribution  to  equation  (58)  from  collisions  between  the 
perturbers  and  the  radiating  atomic  system  is  given  by 


r  U  .  (to) 

a  b  .  a  b 


ijE  £  23  |  / 


a  p  VS  P  a"  p" 


(w-U,,  +  01.-0)  ,,  +  w)  P 
a  b  p  P  p 


+  h  2  E  E  i  / 

b"^b  pM  p  \ 


b  P  I  vs  p  b"  p" 


x  6  (oj  —  o)  +  oo  „  -  o)  +  0)  ,,)  p 
a  b  P  P  p 

+  ^  ^  ^  I  ^ap|vSP|ap^>  -  ^b  p  |v: 

*  P  P'  '  '  N 


(67) 


b  p' 


x  6  (oj  —  o)  +  ov_o),+o))  p. 

a  d  p  p  p' 


We  now  treat  electron  collisions  in  the  binary-collision  approximation,  and  set 

o)  *  o)  which  corresponds  to  the  impact  approximation  in  the  theory  of 
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spectral  line  broadening  . 

It  is  conventional  in  electron-atom  collision  theory  to  introduce  the 
inelastic  scattering  cross  section 
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",  pi  =  -Y~  ^  |  <  ap*]  Vc|a"  ?'•)> 
-ft  V  p"  '  ' 


a  (a  ->a 


and  the  elastic  scattering  amplitude 


6  (a)  -0)  „  -  u)  it  +  to)  (68) 
a  a  p  p 


f  a(P,  P") 


m 


2nfi 


<\  3  P*  I  V°  |  a  P’^  , 


(69) 


where  p  denotes  the  wave  vector  describing  the  electron  scattering  by  the 
.  c 

atomxc  system, and  V  is  the  appropriate  collisional  interaction.  These 
quantities  can  be  defined  without  the  use  of  lowest-order  perturbation  theory 
by  replacing  V  by  the  scattering  matrix  T. 

The  isolated-line  width  produced  by  electron  collisions  is  then  obtained 
in  the  form 


a  b 


£  a  (a-**",  p)  +  V  0  P> 

a'^a  6  a"tb  e 


+  V 

e 


f  (P,  P")  -  f  (P,  P")  2 

a  b 


> 


(70) 


where  Ng  is  the  electron  density  and  Vg  is  the  electron  velocity.  The  density 

matrix  (p)  corresponds  to  the  electron  velocity  distribution.  This  result  is 

1-3 

in  agreement  with  the  work  of  Baranger  .  Note  that  the  total  width  due  to 
collisions  is  the  sum  of  the  rates  for  all  inelastic  transitions  and  a  term 
involving  the  square  of  the  difference  between  the  elastic  scattering  amplitudes. 

We  have  shown  that  the  approximation  of  R(z)  by  /2/*implies  that  the  isolated 
line  width  can  be  obtained  by  adding  up  the  partial  widths  arising  from  auto¬ 
ionization,  spontaneous  and  induced  radiative  transitions,  and  collisions,  as 
expressed  by  equations  (60),  (64),  (66),  and  (70). 
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r 

a  b ,  a  b 


A  +  fSR  +  pi  R  + 

ab,  ab  ab,  ab  ab  ,  ab  a  b,  a  b 


(71) 


Although  the  final  results  have  been  worked  out  explicitly  in  the  approximation 

w  =  -  wb,  it  should  be  emphasized  that  the  general  theory  predicts  that  each 

partial  width  will  have  a  frequency  dependence. 


5.  RELATIVE  IMPORTANCE  OF  VARIOUS  BROADENING  MECHANISMS 

In  order  to  determine  the  relative  importance  of  various  broadening 
mechanisms  for  a  particular  spectral  line,  the  expressions  obtained  for  the 
isolated  line  widths  in  section  IV  must  be  evaluated  using  realistic  atomic 
wavefunctions .  In  this  section,  we  discuss  approximate  evaluations  of  the 
widths  which  reveal  their  dependence  on  the  states  a  and  b  comprising  the  line 
I  and  on  the  physical  properties  of  the  surrounding  particles  and  radiation  field. 

I 

!  The  shifts,  which  are  not  considered  in  this  investigation,  may  be  important, 

l 

i  particularly  in  the  case  of  induced  transitions. 

1 

I 

Even  when  the  atomic  lines  are  assumed  to  be  isolated,  the  upper  and 
lower  levels  a  and  b  usually  consist  of  degenerate  magnetic  sublevels.  The 

!  result  of  superimposing  the  contributions  from  each  allowed  magnetic  component 

I 
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can  be  shown  ’  to  be  the  same  as  the  expression  obtained  from  equation  (45) 

y  2 

after  both  I  D  ,  I  and*X  ,  ,  (uj)  are  summed  over  the  magnetic  sublevels  of 

|  1  a  d  1  a  b ,  a  b 

b  and  averaged  over  the  magnetic  sublevels  of  a.  The  density  matrix  Pfl  must 
be  replaced  by  the  number  density  of  atoms  in  the  state  a,  regardless  of  the 
magnetic  sublevel  quantum  numbers.  In  the  remainder  of  this  section,  all  atomic 
states  will  be  specified  by  giving  only  the  principal  and  angular  momentum  quantum 
numbers  of  the  active  electron. 
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A.  Autoionization 


Consider  the  autoionizaton  process  corresponding  to  the  ejection  of 
an  nil -electron  accompanied  by  the  deexcitation  n  £  n  £  of  one  of  the 


a  a 
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remaining  bound  electrons.  Using  quantum-defect  theory  ,  the  autoionization 

rate  for  large  n  can  be  related  to  the  cross  section  for  the  excitation 

22 

nc£c  na*’a  induced  by  electron  impact.  In  a  previous  paper  ,  we  showed  that 
the  autoionization  rate  due  to  the  dipole  part  of  the  electrostatic  interaction 
can  be  approximated  by 


A  (n  £  ,  n£  -*■ 
a  a  a 


(A) 


■  f  (ft)(5) 

V 

g  (  e  £' 


2<2£  +  1) 

c 


f  (n  £ 
c  c 


n  £  ) 
a  a 


2(2£a  +  1)  2  (2£  +  1) 


ea£) 


(72) 


where  f(n  £  -*n  £  )  is  the  oscillator  strength  and  g  (e  £'  -►  e  £)|  0  is  the 

C  C  3  fl  C  3  £  B 

a 

threshold  value  of  the  partial-wave  Gaunt  factor  for  the  free-free  transition 

of  an  electron  in  the  field  of  the  ion  with  residual  charge  Z.  In  the  non- 

relativistic  approximation  employed  here,  equation  (70)  is  practically  inde- 

-3 

pendent  of  Z.  The  n  dependence  which  is  predicted  is  a  well-known  property. 
The  £-dependence  is  difficult  to  derive,  but  experience  in  evaluating  the  free- 
free  Gaunt  factor  indicates  that  the  autoionlzatlon  rates  decrease  rapidly  with 
increasing  £,  typically  like  £ 

B.  Spontaneous  Radiative  Transitions 

The  properties  of  radiative  transition  rates  are  so  well-known  that 

only  a  brief  discussion  will  be  necessary.  Using  the  asymptotic  behavior  of 

the  hydrogenic  oscillator  strengths  for  large  n,  the  spontaneous  radiative  decay 

23 

rate  due  to  all  downward  transitions  can  be  estimated  by 
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(73) 


Ar(n) 


1.6  x  10 


10 


-1 

sec 


It  should  be  kept  in  mind  that  equation  (73)  describes  only  An^  o 

4 

transitions,  for  which  the  radiative  decay  rates  increase  like  Z  with  increasing 
-9/5 

Z.  The  n  dependence  indicates  that  the  dominant  radiative  decay  mode  for  a 
doubly-excited  state  with  the  outer-electron  in  a  high-n  level  will  be  the 
radiative  deexcitation  of  the  inner-electron. 

The  total  spontaneous  radiative  width  is  estimated  by 


„SR 


a  b  *  a  b 


=  1.6  x  10 


10  4 


(n 


-9/2 


-9/2, 

%  > 


(74) 


with  the  understanding  that  n^  and  n^  refer  to  the  lowest  values  for  which  a 
dipole  transition  can  occur. 

In  medium- and high-Z  ions,  the  most  prominent  emission  lines  have  radiative 
decay  rates  which  are  at  least  comparable  to  any  autoionization  rates.  Con¬ 
sequently,  equation  (74)  can  be  used  as  an  estimate  for  the  total  natural 
width.  It  should  be  emphasized  that  prominent  absorption  lines  can  arise  from 
transitions  to  autoionizing  levels  whose  spontaneous  radiative  decay  rates  are 
relatively  small. 


C.  Induced  Radiative  Transitions 

The  result  obtained  after  averaging  equation  (66)  over  the 
quantum  numbers  can  be  written  in  the  single-particle  approximation 


magnetic 

as 
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We  note  that  the  electron  Impact  width  associated  with  inelastic  transitions 

can  be  obtained  in  the  Bethe  approximation  from  equation  (75)  simply  by 

Ilf  2  -  f  (V  )  - 

replacing  -  x  J  I^w.A)  d/V  by  e  Ng  j  “■  . 

e  \{3 

2  24 

The  application  of  the  r  sum  rules  given  by  Bethe  and  Salpeter  is 

prevented  by  the  photon-frequency  dependence  of  the  radiation  field  intensity 

1^  (u,ft).  If  we  define  a  total  intensity  by 


I  (u>) 


?/ 


dJV  <“"V 


(76) 


and  replace  I(w)  by  a  suitable  frequency-independent  average  intensity  I,  the 
summations  over  the  complete  sets  of  states  n  „l  „  and  n.  „JL  „  can  be  carried 

&  &  D  D 

out,  and  the  absorption  and  stimulated  emission  contributions  can  be  combined  to 
give  the  result 
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The  expectation  values  of  r  will  now  be  approximated  by  using  the  hydrogenic 
formula 
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/  I  2 1  ^ 

n£.  |  r  | 


2Z 
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(78) 


2n  2 

=  7"  ^ 


where  the  last  version  corresponds  to  an  it-independent  approximation. 

Let  us  consider  a  thermal  equilibrium  distribution  of  radiation  character¬ 
ized  by  a  radiation  temperature  T^.  We  assume  that  this  radiation  is  incident 
from  one  side  only,  so  that  the  integration  over  0  ^  reduces  to  multiplication 

by  2  7t.  We  choose  the  frequency-independent  intensity 

.3 


I  = 


(kTR) 

2w2cV 


(79) 


which  is  about  2/3  the  maximum  value  of  the  Planck  distribution  function 

corresponding  to  TD.  The  estimated  width  due  to  induced  radiative  transitions 
R 


is  then  given  by 
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(80) 


2.67  x  10‘ 


sec 


-1 


SR  I  R 

In  contrast  to  equation  (74)  for  r  ,  .  ,  T  .  .  increases  with  increasing 

ab, ab  ab, ab 

principal  quantum  number  and  decreases  with  increasing  residual  charge.  The 

(k  T_)  dependence  on  the  radiation  temperature  is  also  an  important  property. 
R 
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Unless  the  Induced  radiative  transitions  have  frequencies  within  the  region  of 
the  maximum  in  I(ui),  equation  (80)  is  probably  an  overestimate  of  the  actual 
width. 


D.  Electron  Collisions 

25 

It  has  been  pointed  out  by  Griem  that  the  Bethe  approximation  for 
the  electron-impact  width  due  to  inelastic  transitions  can  be  modified  to  take 
into  account  some  of  the  effects  of  elastic  scattering.  This  is  done  by 
beginning  the  integral  over  the  Maxwellian  distribution  at  zero  velocity  rather 
than  starting  it  at  the  excitation  threshold.  The  result  can  be  written  as 


E 

n"  i"in  l 
a  a  1  a  a 
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V  V+Vb 
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where  g  now  denotes  the  thermal  average  of  the  free-free  Gaunt  factor. 

If  the  Gaunt  factors  in  equation  (81)  are  now  replaced  by  suitable  average 
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values,  the  r  sum  rules  can  be  employed  to  obtain  the  result 
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(82) 


b  ~b 


where  AE^  and  AE^  should  be  chosen  to  be  the  smallest  threshold  energies.  We 
now  use  equation  (78)  and  set  the  average  Gaunt  factors  equal  to  unity,  which 
is  a  good  approximation  for  An=0  transitions  in  multiply-charged  ions.  The 
electron  impact  width  is  then  estimated  by 


(83) 


As  expected,  the  electron  impact  width  has  the  same  n-  and  Z-  dependences 

as  the  width  due  to  induced  radiative  transitions.  The  Jl-dependences  would 

also  be  the  same  if  they  had  been  retained  in  the  approximation  for  the 

2 

expectation  values  of  r  .  The  condition  for  thermal  radiation  broadening  to  be 
dominant  over  electron  impact  broadening  is  found  from  equations  (80)  and  (83) 
to  be 


(84) 
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